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ABSTRACT 

Context. Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary 
structure commonly formed before and during the formation of cores. 

Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. 
Methods. We use Herschel observati ons of the Serpens Core to compute temperature and column density maps of the region. Among 
the simulations of [Dale et al.| ( |2012l ), we select the early stages of their Run I, before stellar feedback is initiated, with similar total 
mass and physical size as the Serpens Core. We derive temperature and column density maps also from the simulations. The observed 
distribution of column densities of the filaments has been analysed first including and then masking the cores. The same analysis has 
been performed on the simulations as well. 

Results. A radial network of filaments has been detected in the Serpens Core. The analysed simulation shows a striking morphological 
resemblance to the observed structures. The column density distribution of simulated filaments without cores shows only a log-normal 
distribution, while the observed filaments show a power-law tail. The power-law tail becomes evident in the simulation if one focuses 
just on the column density distribution of the cores. In contrast, the observed cores show a fiat distribution. 

Conclusions. Even though the simulated and observed filaments are subjectively similar-looking, we find that they behave in very dif¬ 
ferent ways. The simulated filaments are turbulence-dominated regions, the observed filaments are instead self-gravitating structures 
that will probably fragment into cores. 

Key words. Stars: formation; ISM: structure; ISM: Evolution, ISM: individual (Serpens Main, LDN 583, Ser G3-G6) 


1. Introduction 


Eilamentary molecular structures are commonplace in giant 
molecular clouds (GMCs). Large-scale filaments have been 
detected in extinction maps of low-mass star forming regions 
(e.g. |Cambresy||1999[ [Kainulainen et al.||2Qll| ) or in the sub- 
/millimeter wavelength range in high-mass star forming regions 
(e.g. [Johnstone & Bally||1999t [Gutermuth et al.|[2Q08 ). In the 
last years observations have been carried out using the Herschel 
telescope in nearby star-forming regions, as well as in regions 
where star formation is not currently active (e.g. the Polaris 
Elare region, [Andre et~ar]|2010| ). A characteristic filamentary 
structure consists of a main filament with some sub-filaments 
converging on it. This is the case in low-mass star-forming re¬ 
gions (e.g. the B213 filament in Taurus, or the Corona Australis 
region; [Palmeirim et al.[[?013t [Sicilia-Aguilar et al.[[MT3] ), as 


^ Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


well as in high-mass star-forming regions (e.g. the DR21 ridge 
and filaments in Cygnus; [Hennemann et al.|2Q12 1. 

Eilaments are thought to be the consequence of thermal, dynam¬ 
ical, gravitational or magnetic instabilities. These can be caused 
by collisions of large-scale warm flows that are themselves gen¬ 
erated by the turbulent velocity fields that are also ubiquitous 
in GMCs. The filaments originated by gravitational instabili- 
ties can be either produced i n dense two-dimensional structures 
( [Burkert & Hartmann|2004[) or resulting from the expansion of 
a feedback-driven bubble ( [Krause et al.|2013[ ), or from the col¬ 
lision of two GMCs. An extensive review on the possible mech¬ 
anisms of filament formation is presented in [Andre et al. [ ( [20Id) 
and [Molinari et al.[ ( [2014[ ). 

Numerical simulations of star formation on GMC-scales have 
in recent years become sophisticated enough that detailed com¬ 
parisons with observations, beyond simple statistics such as star 
formation efficiencies and stellar mass functions, are required to 
move forward. Attempting to model a given object or system 
in detail is fraught with peril, since the initial conditions from 
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which to begin the model are of course not known, and even the 
present-day state of the system is unlikely to be understood in 
enough detail. However, many simulations appear to be able to 
reproduce general features of star-f orming clouds, such as fil - 
amentary structures in the gas (e.g. |Moeckel & Burkert||2Q14| ). 
To determine if this resemblance is more than merely subjective 
requires state-of-the-art observations of candidate systems and 
detailed comparisons with suitable simulations. 

Recently, a comparison between GMC simulations and obser¬ 
vations has been presented by [Smith etH] ( |2014| ) in order to 
investigate the general observational result that filaments appear 
to have a surprising constant width of O.lpc regardless of the 
structure or properties of their host cloud ( [Arzoumanian et ak] 
[2011[ Andre et al.|2014[ ). The authors performed hydrodynami- 
cal simulations of molecular clouds with turbulent velocity that 
is responsible for complex filamentary fields. They found that 
the widths of the filaments depend on the data range fitted by 
a Gaussian function, with an average full width half maximum 
(FWHM) of 0.3 pc. This value is in agreement with theoretical 
predictions for accreting filaments ( [Smith et al.|2014[ ). 

A different approach, presented here in this paper, consists in the 
comparison between a specific GMC numerical simulation from 
[Dale et alT]p012[ ) and a particular real object, namely a region of 
the Serpens molecular cloud. 

The Serpens molecular cloud represents an ideal nearby labora¬ 
tory to study on-going low-mass star formation. The distance 
to the Serpens Core of 415+5 pc has been recently revised with 
astrom etric observations obtained with the VLBA ( [Dzib et ak] 
[2010[ ). The authors also estimated a distance for the entire Ser¬ 
pens molecular cloud of 415+25 pc. 

The region is located north of the Galactic plane and part of the 
Aquila Rift. [Konyves et al.[ ( |MT0| ) and [Bontemps et al^ ( [2010[ ) 
presented the analysis of the Herschel observations of the Aquila 
Rift molecular complex and Serpens South, which are both in 
the southern part of the Serpens molecular complex. Most of 
the protostars resolved in the Herschel observations were newly 
discovered. In the Aquila rift the spatial distribution of protostars 
indicates three star-forming regions, the richest of which is the 
W40/Sh2-64 HII region. 

The new Herschel maps of Serpens presented in this paper are 
dominated by the so-called ‘Serpens Core Cluster’ and the com¬ 
plex around a group of four T T auri stars, called G3 to G6 (see 
Figure [T]). In their Spitzer study Harvey et al.[ ( [2006| ) called the 
Core Cluster ‘Cluster A and the southern cluster around Ser G3- 
G6 ‘Cluster B’. Located further south, t he prominent Herbig Be 
star VV Ser ( [Harvey & Dunharn|[2009[ ) is visible in the 70 jim 
Herschel map as a point source. As highlighted in the review 
of [Eiroa et al.[ ( [200i^ , a peculiarity of the Serpens Core is that 
all the phenomena related to active star formation have been 
discovered, from still in-falling gas to pre-protostellar gaseous 
and dusty condensations, as well as already-formed Class 0 to 
Class II young stars. A similar scenario is found also in other 
young star-forming regions as, for example, the Coronet cluster 
([Sicilia- Aguilar et al.[[2011[), or IC1396A ([Sicilia-Aguilar et al.[ 

[2014[[2015[ ). _ _ 

In the millimeter regime [Enoch et al.] ( [2007[ ) detected sources in 
Cluster A, but also in Cluster B and in the bright filament north¬ 
east of it. Several sub-millimeter sources have been detected in 
the Core Cluster. They are associated with a south-eastern and a 
north-western clump ( [Davis et al.[1999t[Casali et al.[1993[ ). Con¬ 
tinuum sub-millimeter observations carried out with the SCUBA 
bolometer array camera on the James Clerk Maxwell Telescope 
(JCMT) at 450 and 850yum revealed a prominent diffuse fila¬ 
mentary structure in the south and east direction from the south¬ 


eastern cluster of sub-millimeter sources. The outflows driven by 
these objects have been locally constrained close to the sources 
by C O observations ( [White k al.[[1995t [Duarte-Cabral et ak] 
[2010[ ). Observations of molecular gas obtained with the IRAM 
30m telescope revealed two velocity components between the 
two sub-clusters in the Serpens Core ( [Duarte-Cabral et al.[2010[ ). 
In their paper the authors proposed that after an initial trig¬ 
gered collapse, a cloud-cloud collision between the northern and 
southern clouds is responsible for the observed characteristics 
of the two sub-clusters. Observations of A 2 H^, HCO^, and 
HCN toward the Serpens Main molecular cloud have been re¬ 
cently obtained using CARMA ( [Lee et al^[2014[ ). The authors 


find that while the north-west (NW) sub-cluster in the Serpens 
Core shows a relatively uniform velocity, the south-east (SE) 
sub-cluster shows a more complicated velocity pattern. In the 
SE sub-cluster they also identified six filaments: those located 
in the northeast of the SE sub-cluster, have larger velocity gra¬ 
dients, smaller masses, and nearly critical mass-per-unit-length 
ratios, while the opposite properties have been found in the fila¬ 
ments in the southwest of the SE sub-cluster. 

In this paper we present the continuum Herschel observations 
of the Serpens Core and compare them to a numerical simula¬ 
tion performed by [Dale et al.[ ( [2Q12| ). We organize the paper as 
follows: in Section [2| we summarize the Herschel observations 
and the data reduction. In Sectionwe present the analysis car¬ 
ried out including the derivation of the temperature and the col¬ 
umn density maps, while in Section [4] the simulations used in 
our study are described. In Section [^we compare the proper¬ 
ties of the observed filaments with the filaments formed in the 
simulations. A summary and conclusions are given in Section 

2. Observations & data reduction 

The Herschel data of Serpens were taken from the Herschel Sci¬ 
ence Archive (HSA) and are part of the Herschel Gould Belt sur¬ 
vey ( [Andre et al.|2Q10[ ). The observations (obsIDs: 1342206676 
and 1342206695) were made on October 16th, 2010 and cover 
an area of 2.5 deg x2.5 deg. The Serpens cloud was simulta¬ 
neously imaged with Herschel ( [Pilbratt et al. 2010[) in fiv e wave¬ 
lengths, using the two cam eras PACS ([Poglitsch et ar[[2010[ ) at 
70 and 160yum and SPIRE ( [Griffin et al.|2Q10[ ) at 250, 350, and 
500 yum. Two orthogonal scan maps were obtained by mapping 
in the parallel fast scan mode with a speed of 60'7s. We reduced 
the data using HIPE vl0.3.0 for the calibration of the Level 0 
PACS and SPIRE data. Then the Level 1 data of both instru¬ 
ments were taken t o produce the fi nal maps with the Scanamor- 
phos package v23 ( [Roussel|2013[ ), using the jparallel option for 
the reduction of observations with PACS and SPIRE in parallel 
mode and the Igalactic option to preserve brightness gradients 
over the field. To properly take care of glitches and artefacts that 
occurred in the PACS maps, we canceled the HIPE glitch-mask 
at the conversion of the Level 1 data to the Scanamorphos for¬ 
mat. Additionally, we switched ihtijumps_j>acs option on for the 
production of the final maps in order to first detect and then re¬ 
move the jumps in the PACS mosaics. Those discontinuities can 
be caused either by glitches or by electronic instabilities. The 
pixel-sizes for the five maps at 70, 160, 250, 350, and 500yum 
were chosen as 3.2", 4.5", 6", 8", and 11.5", respectively. The 
recent version of Scanamorphos, in particular, takes into account 
the PACS distortion fiat-fields. 

The final mosaics have been calibrated applying the Pl anck off- 
sets. We adopt the rigorous approach presented by e.g. [Bernard[ 
[et al.[ p010[ ) which consists in combining the Planck data with 
the IRAS data of the Herschel fields and extrapolating the ex- 
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Table 1. Planck offsets of the Hersch el mosaics of the Ser pens Core 
obtained using the method presented in [Bernard et al.| ( [T010| ). 


Instrument 

Filter 

[jum] 

Offset 

[MJy/sr] 

PACS 

70 

12.2 


160 

111.6 

SPIRE 

250 

74.9 


350 

37.6 


500 

15.8 


pected total fluxes to the Herschel filter bands. The difference 
between the expected and the observed fluxes provides the off¬ 
sets to correct the Herschel mosaics. The Planck offsets obtained 
at each filter are presented in Table These values are added to 
the Level 2 mosaics obtained with Scanamorphos. 


3. Analysis 

In this section we present the analysis carried out on the Her¬ 
schel o bserv ations of the Serpens molecular cloud. In particular, 
Section [3?T] presents the morphological description of the “Ser¬ 
pens C ore” (Cluster A) and the cluster “Ser G3-G6” (Cluster B). 
Section [L2| presents temperature and column density maps of the 
region. 

3.1. Molecular cloud morphology 

Figure[2shows the composite of Herschel/PACS 70yum, 160 pm, 
and SPIRE 350yum. All the cloud structures become brighter at 
longer wavelengths. The entire Herschel maps in the eastern 
part extend up to the far-infrared bright region LDN 583. This 
region appears to be almost empty, with only a cloud structure, 
prominent at the SPIRE wavelengths, with a point source, in the 
South-West direction from LDN 583. 

The Serpens star-forming region shows a highly filamentary 
morphology. The Serpens Core Cluster seems to be connected 
by several filaments to the underlying molecular cloud. These 
filaments have a radial structure similar to spokes of a wheel 
with the cluster forming the hub. This morphology is a common 


feature of star-forming regions (e.g. 

Myers 

2009 [ [Hennemann[ 

et al.[[2012). [Schneider et al.[([2012 

) note t 

lat the junction of 


filaments seems to be intimately connected with the formation 


of star clusters, as was also found by [Dale & Bonnell| ( |2011| ) in 
numerical simulations, suggesting that gas flows along the fil¬ 
aments and accumulates at the hubs. The filaments are very 
structured showing kinks and condensations. In particular, the 
westernmost filament harbors a bright compact source. With 
respect to the colors, the eastern filaments show a remarkably 
bright emission at 70yum at their tips. 

The Ser G3-G6 cluster seems to be aligned with a filament and 
associated with several bright sources in its vicinity. The fila¬ 
ments or cloud edges seen in this region are highly structured 
and reminiscent of turbulence. 


3.2. Temperature and column density of the clouds 

The cloud temperature has been obtained using two different ap¬ 
proaches: the first approach uses the ratio of the fluxes at 10 pm 
and 160/im, while the second is based on the fit of the SED 


fluxes between 70 and 500 prn^ In both cases we used the Her¬ 
schel maps after applying the Planck offsets. While the first ap¬ 
proach conserves the angular resolution of the 160 pm image (as 
done in [Preibisch et al.||2012HSicilia- Aguilar et al.| 2014|), the 


second uses the PACS and SPIRE images all convo^ 


ved to the 


500 yum resolution (as done in [Roccatagliata et al^|2013| ). The 
fit of the fluxes is obtained by using a black-body leaving as free 
parameters the temperature T and the surface density E [g/cm^], 
with the same procedure presented in [Roccatagliata et aL ( [2013| ). 
The dust properties are fixed by interpolating at the Herschel 
bands the values of the dust mass absorption coefficient Ky from 
the model of [Qssenkopf & Henning[ ( p^94| ) with a value of 1.9, 
which adequately describes a dense molecular cloud. As already 
found and commented by 


fitting leads to slightly co; 


Roccatagliata et al.[ ( [2013| ), the SED 
he temper- 


ilder temperatures. Since t 

ature T and the surface density E [g/cm^] are two independent 
parameters, the best fit parameters have been obtained minimiz¬ 
ing;^^. Eor the further analysis presented in the paper we use the 
temperature and surface density maps obtained with this SED- 
fitting method. 

The column density is computed using the relation 


Ah = 2Ah, = 


2SR 
mn • PYi2 


( 1 ) 


where S is the surface density, mn is the hydrogen mass (i.e. 
L6736e-24 g) and p^^^ is the mean molecular weight (i.e. 2.8). 
Multiplying by the gas-to-dust mass ratio R (assumed to be 100), 
we obtain the total column density. 

In the entire region covered by Herschel observations we find 
that the average temperature of most of the nebula is about 24 K. 
It is interesting to notice that the coldest part of the cloud, of 
about 19 K, corresponds to the filamentary material converging 
into the central core. 

The total mass for each region considered in our study is com¬ 
puted summing up the column densities of all the pixels. Us¬ 
ing a distance of the Serpens molecular cloud of 415+25 pc, 
the total mass (of dust and gas) over a region of 2.0° x 2.2° is 
3213 ^ 3 ^^ Mq. The errors have been computed considering the 
error of 25 pc on the distance ( [Dzib et ar][201Q[ ), which repre¬ 
sents the dominant error in this estimate. 

The total mass of the Ser G3-G6 clusterE]and the LDN 583 re¬ 
gion His 98^12 Mo and 67 Mq, respectively. The temperatures 
rangM^etween 11 K and 16 K in the case of the Ser G3-G6 clus¬ 
ter, while in the case of the LDN 583 region they range between 
14 K and 18 K. 


A zoom of the temperature and column density map of the Ser¬ 
pens Core is shown in Eigure The total mass present in the 
Serpens Core is 494 Mq, where the errors are computed from 
the error in the distance. The region considered for the total mass 
of the Serpens Core is highlighted with a large box in Eigure 
and is 42'x42', which corresponds to about 5 pc x 5 pc. 

Using the extinction thresholds in the study of [Lada et 
\20\2\ , we computed the following total masses in the re¬ 
gion of 5 pc X 5 pc of the Serpens Core: for Ak > 0.1 mag 


^ We highlight that this method can have some caveats in using also 
the 10 pm map, since the emission might trace a hotter dust population 
compared to the colder dust mapped at longer wavelengths. 

^ computed in a box of 7.3' x 13.4' in size, center (crj 2 ooo. <ij 20 oo) = 
(18^ 28”^ 58% +00° 37'55"), and orientation 37° in North-East direc¬ 
tion. 

^ computed in a box of 21.7' x 9.3' in size, center (crj 2 ooo. ^j 2 ooo) = 
(18^ 33"^ 37% +00° 30'07"), and orientation 40° in North-East direc¬ 
tion. 
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Fig. 1. Composite HerscheljVACS lOjjm (blue), 160 jim (green) and HerscheljSVWE 350yum (red) mosaics in logarithmic scale. The ranges 
are between 0 and 0.3 Jy/arcsec^, 0 and 0.3 Jy/arcsec^, and 0 and 0.05 Jy/arcsec^ for the blue, the green, and the red image, respectively. The 
center of the figure corresponds to the position (crj 2 ooo, <^j 2 ooo) = (18^ 30"^ 49% +00° 39' 24") and the size of the region is 1.7° x 1.9°. 


{Ay ^ 0.9 mag), the total cloud mass is and for 

Ak > 0.8 mag {Ay ^ 7.3 mag) IOO^JqMq. The fraction of 
dense material present in the Serpens Core (using as threshold 
Ak > 0.8 mag) is ~ 20%. The star formation rate o f the Ser¬ 
pens Core has been computed by [Harvey et al.| ( [2007| ). Rescal¬ 
ing their value with the updated distance of 415 pc of the Serpens 
Core, we obtain 28 Mo/Myr. These values for the Serpens Core 
agree well with the relation derived for other low-mass molec- 
ular cl ouds by [Lada et aT] ( |2012| |2010| ) and [Heiderman et al 
(|2010|). This rel ation has been further debated by |BurkerF& 


Hartmann] pO 13 j ) , who do not find a particular threshold for star 


the Serpens molecular cloud of 311 pc ( jde Lara et al.||199T] ), 
[White et aT] ( |1995t estimated a lower mass limit of 1450 M© 
from the C^^O (J=2-l) data and 400 M© from the C^^O data. We 
scaled their masses to the new distance of 415 pc and computed 
the masses of the same regions from our data. We found a total 
mass about 7 times lower than the gas mass estimate from the 
C^^O. A possible reason for this discrepancy might be a different 
value of the gas-to-dust ratio which is always assumed to be 100 
in our calculations. 

Testi & Sargent[ ( p^98] ) computed the total mass from the 3 mm 


formation but instead a strong continuous increase of the star 
formation efficiency with density. 


3.2.1. Comparison with previous mass determinations: dust 
vs. gas masses 

We compare the total mass computed for the Core region with 
previous dust and gas mass estimates. Assuming a distance of 


continuum emission in the inner part of the Serpens Core within 
5'x5'. Out of the 4 cores detected by [Testi & Sargent|([1998| ), 
only the mass of 3 of them is reported in [Testi efaL p000[ ). We 
scaled their masses to the 415 pc distance and we computed the 
masses using our column density map on the same positions of 
their cores. Only Core D, which corresponds to the position of 
S68N, has the same mass using both approaches. Cores B and 
C are both three times less massive than the masses derived at 
3 mm. [Lee et al^ ( [2014[ ) presented the properties of 18 cores 
detected in the Serpens Core in the 3 mm continuum data ob- 
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tained with CARMA. The Core B and Core C in |Testi & Sargent 


(1998) correspond to the cores Sll and S7, respectively, in 


Lee 


et al.| ( 20T4) . They computed a mass of about bM© for Sll and 


11 M© for S7. We highlight that the mass of Sll from |Lee et al.| 
( 2014| ) corresponds to the mass we obtained. The mass of S7 
differs only marginally, while other cores are about two times 


Testi et al. 

(I 2 OOOI). 

Duarte-Ca 

3ral et al. 


masses of the clumps identihed in the Serpens Core. Th e clumps 
ontheC^^Q J =1-0 channel maps have been derived by 


Duarte- 


Cabral et al. (2010) using the 2DCLUMPFIND algorithm, as¬ 


suming a constant temperature of 10 K, a mean molecular weight 
of 2.33, and a C^^O fractional abundance of 4.7x10“^ with re¬ 
spect to H 2 . They computed also the values of the virial mass 
of the clumps from the observed velocity dispersion. The ra¬ 
tio between the virial mass (Myiriai) and the gas clumps mass 
(^gas) gives the information whether the structure is gravitation¬ 
ally bound (Mviriai/4^gas:^l) or unbound (Myiriai/^^gas >!)• Ac¬ 
cording to this dehnition, cores D, E, F, and G are unbound. 
Ah the dust masses computed from our Herschel maps in corre¬ 
spondence to the clump positions agree well with the estimates 
based on the sub-millimeter and millimeter fluxes. We highlight 
that the total mass that we compute for the dust mass assumes 
that the gas-to-dust ratio is 100. However the adopted canonical 
value has an intrinsic uncertainty of at least a factor of 2, espe¬ 
cially locally. This can explain the small differences found in the 
presented comparison. 


4. Star formation simulations. 


Dale et aL] ( |2012| ) simulated the evolution of a group of model 


GMCs distributed in (mass, radius, velocity dispersion) param¬ 
eter space so as to mimic the properties of Milky Way clouds 
inferred by |Heyer et al.| ( [200^ . The clouds range in mass from 
lO'^-lO^ M© and in initial radius from 5-90 pc. The initial condi¬ 
tion of the simulations is an idealised spherical cloud with mild 
Gaussian density prohles giving a degree of intrinsic central con¬ 
centration, seeded with a supersonic turbulent velocity held. The 
velocity held initially has the power spectrum of Burgers turbu¬ 
lence with P{k) oc and wave numbers between 4 and 128 are 
initially populated. The normalisation of the velocity held is ad¬ 
justed to give ah the clouds an initial virial ratio of 0.7, so that 
they are ah formally bound. This results in a range of velocities 
of 2-10 km s“^. 


The density held responds rapidly to the velocity held and the 
clouds develop complex structures of sheets and hlaments from 
the interactions of turbulent shocks. Despite being formed by 
turbulence, many of these structures are persistent, although 
they are often advected through the clouds by the hows from 
which they formed. On timescales comparable to the cloud 
freefah times (~ 1 Myr), sufficient mass accumulates in some 
dense structures for subregions of the clouds to become self- 
gravitating and to begin forming stars. 

The purpose of the simulations was to examine how the effects 
of massive stellar feedback in the form of photoionisation and 
winds vary with the properties of the host clouds. However, 
these calculations have other uses. In particular, one can search 
through the parameter space (including the time dimension) to 
hnd a simulation, or a part of a simulation, which is similar to 
a given real system. Serendipitous similarities are likely to be 
more trustworthy than targeted simulations because no assump¬ 
tions have been made to achieve a desired result and one can 


be conhdent that the result is not there merely because it was 
inserted from the beginning. 

Amongst the calculations of [Dale et al.| ( |2Q12'| ), we hnd that the 
innermost 5pc of their Run I, at a time that we refer to as t^c 
(0.25 Myr after the initiation of star formation but before stellar 
feedback is initiated), bear a striking (subjective) morphological 
resemblance to the Serpens Core. The gross physical character¬ 
istics such as mean column and volume densities, and temper¬ 
atures are also numerically sufficiently similar to make a more 
detailed comparison worthwhile. We therefore treat this subre¬ 
gion of the Run I calculation as a model of the Serpens Core, 
in the same way that the real Serpens Core is a subregion of its 
parent molecular cloud. 

A 5 X 5 X 5pc box (shown in projection in Figure 3) near the sim¬ 
ulation center of mass at this epoch exhibits a pc-scale core con¬ 
taining a small cluster of three stars with a total mass of 14 M©, 
with several dense hlaments radiating from it, one of which joins 
with a smaller core approximately one pc distant in projection 
(and in three dimensions, in fact). The masses of the larger and 
smaller cores are ~ 100 M© (Box 1 in Figure and ~ 30 M© 
(Box 2 in Figurej^, and total gas mass column and volume den¬ 
sities are Ig cm“^ and ~ lO'^-lO^cm"^, respectively. There are 
also several other hlamentary structures connected with the two 
cores. 

The hlaments in the simulation were originally created by col¬ 
liding hows in the supersonic turbulent velocity held imposed as 
part of the simulation initial conditions. The clouds in the simu¬ 
lations of [Dale et al.H2012| are ah initially mildly centrally con¬ 
densed and gravitationally bound, so that the densest structures 
tend to form at the bottom of the clouds’ gravitational wells near 
their barycenters. However, owing to the different mean densi¬ 
ties and turbulent velocities of the clouds, the properties of the 
hlaments vary between simulations. Run I is the simulation with 
the smallest mean turbulent velocity and a fairly low mean den¬ 
sity, so the densities of the hlaments are also comparatively low. 
They form a hub structure and material hows along them towards 
the common hub, where it accumulates and forms a dense cluster 
of stars. The timestep, tsc, chosen for this comparison occurs at 
the very beginning of the formation of the cluster, when there are 
only a few stars present. We also use a second snapshot 0.4 Myr 
further forward in time to examine how the characteristics of the 
central region of the simulation vary on timescales of the order 
of the freefah time, in particular the column density distribution. 
In contrast to what occurs in other simulations from IDale et ah] 
( 2012| ), the hlaments in Run I do not become gravitationally un¬ 
stable and do not fragment to form stars themselves because 
the freefah timescale for fragmentation to occur within them 
is longer than the timescale on which hows along the hlaments 
transport material into the hub where they join. They are long- 
lived structures whose function over ~Myr timescales is largely 
to deliver fresh gas to the forming cluster located near the cloud’s 
center of mass. 


4.1. Analysis of the Simulations 

We identify in Figure three hlamentary structures of the Run I 
simulation at the selected time tsc, one of which physically con¬ 
nects the two cores. The masses of the hlaments are in the range 
20-40 M©, comparable to the mass of the small core and a few 
times less than that of the larger core (see Table [^. The how 
velocities of the gas are of the order of 1 km s“^. The hlaments 
do not fragment and most of the gas they contain is accreted by 
the cores (and the clusters they eventually form). It is simplistic, 
however, to describe the gas as only howing along the hlaments. 
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At early times when the cloud velocity field is still dominated by 
the original imposed turbulence, the filaments are being advected 
by larger-scale fiows in directions which are not necessarily par¬ 
allel or perpendicular to their long axes. The velocity field in the 
neighbourhood of the two cores is extremely complex and, at 
this early stage in the simulation, has not yet settled down into 
the simpler configuration of several filaments feeding a common 
hub. 

The larger and smaller cores begin to fragment into stars and 
merge into a single object after IMyr and, after this occurs, 
the cluster continues to grow by accretion in a morphologically 
simpler filament-hub configuration. From a control simulation 
without feedback, it is clear that this morphology would be sta¬ 
ble for a protracted period of time, even while the central pro¬ 
tocluster continues to grow and form stars. It rapidly destroys 
both the core and the filament network, leaving an eroded pil¬ 
lar and a fully-exposed cluster. We analyse the raw temperature 
and column density maps extracted directly from the hydrody- 
namical simulations without performing any radiative transfer 
calculations. We rescale the pixel scale and convolve the surface 
density maps of the simulat ions with the Point S pread Function 
(PSF) of SPIRE at 500yum ( Gordon et al.||2008|). W e apply the 
FluxCompensator (as done in Koepferl et al.|2Q15] ) to perform 
this task|3 We use the same conversion factor as in the observa¬ 
tions to convert projected gas mass surface densities to column 
densities as in Equation 


5. Discussion: observed fiiaments versus 
simuiated fiiaments. 

In order to compare the properties of the filaments in the obser¬ 
vations of the Serpens Core to those formed in the simulation, 
we carry out in parallel the same analysis on the observations 
and simulations. The pixel scale of both, observations and sim¬ 
ulations, is about 0.02 pc/px. Both are convolved to the spatial 
resolution of SPIRE at 500yum. 

In our analysis, we refer to compact enhancements in column 
density as ‘filaments’ if they are long and linear, and as ‘cores’ 
otherwise. Using this definition in the observations, the cores 
always contain point source counterparts at 70yum. The analy¬ 
sis of the observed filaments is carried out by defining by eye 
a box around each filament. The positions of the 9 filaments 
analysed are highlighted in Figure The filaments numbered 
1 to 8 are converging into the Serpens Core, while the ninth is 
placed between the Serpens Core and the Ser G3-G6 region. In 
the appendix we zoom in on each filament, showing the Her- 
schel, the column density, and the temperature maps. In order 
to define the boxes of the cores present in filaments 1 and 7, we 
checked the increase in column density and the emission at 70 
and 160 yum, where the point sources are resolved. The boxes 
which include the cores are much larger than the point sources 
resolved at 70 yum: this conservative approach allows to exclude 
any contamination by the cores in the analysis of the filamentary 
structures alone. 

Five boxes (highlighted in Figure have been drawn on the 
most prominent structures in the simulations. These include the 

^ The FluxCompensator was designed to produce realistic synthetic 
observations f rom radiative transfer modelling (e.g. Hyperion by |Ro-| 
|bitaille|[2011| ). By modelling the effects of convolution with arbi¬ 
trary PSFs, spectral transmission curves, finite pixel resolution, noise 
and reddening, hydrodynamical simulations - through radiative trans¬ 
fer modelling - are directly comparable with real observations ( www. 
inpia-hd.inpg.de/~koepferl/). The tool will be publicly available 
in the future. 
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main central core where a small cluster of four stars has already 
formed, the secondary core, a filament which bridges the gap 
between the main and secondary core, and two other prominent 
filaments. 

In order to understand the nature of the network of filaments 
present in the Serpens Core, we construct the histograms of the 
column densities of all the filaments combined together in or¬ 
der to minimize the error in the histogram (which corresponds 
to V^). In each histogram we report on the upper x-axis the vi¬ 
sual extinction corresponding to the column density, computed 
adopting the canonical relation of, e.g., [Bohlin et al.| ( |1978| ): 
Ay = 1 mag corresponds to Ah ~ 2x10^^ cm“^. 

We also compute the column density pr obability density functio n 
(PDF) following the definition given in Schneider et al. ( 2015| |~' 


5.1. Mass and Temperature of the filaments 


We first compute the total mass and the temperature range in 
each observed region. The results are reported in Table For 
comparison, we also show the characteristics of filament number 
9, which is not related to the network of filaments of the Serpens 
Core. 

The masses range between 2 and 33 M©. Summing up the 
masses of the network of filaments converging into the Serpens 
Core a total mass of gas and dust of 95 M© is computed with¬ 
out considering the massive filament (Number 9 in Fig.[^. This 
represents about 16% of the total mass derived from the col¬ 
umn density map in the Serpens Core region from the Herschel 
wavelengths. The most massive filaments have temperatures a 
few K colder than the less massive ones. Dividing these masses 
by the lengths of each filament, we also computed the “mass 
per length” (Ml). The comparison of Ml with the critical val¬ 
ues of mass length (ML^crit) computed for an isothermal, self- 
gravitating cylinder without magnetic support, allows to under¬ 
stand if the structure is thermally supercritical and possibly in 
gravitational contraction. Following |Qstriker| ( |1964| ) we com¬ 
pute ML,crit as: 

ML,cnt = 2 c]IG = 16.7 M^pc-^ (2) 


where the temperature T is the average temperature of each fil¬ 
ament. Only the filament number 1 has a Ml larger than Mi^crit 
and might gravitationally contract. All the other filaments have 
Ml lower than M^crit^ 

Some of the promine nt filaments analy sed in our work have been 
recently resolved by Lee et al. ( 2014 ) with A2H^(l-0) observa¬ 
tions carried out with CARMA. All their filaments are found to 
be associated with the SE sub-cluster in the Serpens Core, while 
our Herschel coverage allowed us to detect 4 filaments (numbers 
5, 6, 7, and 8) also close to the NW sub-cluster. Their filaments 
ESI, FS2 and FS3 correspond to our filament number 1, while 
their structures FCl and FNl correspond to numbers 2 and 4, 
respectively. Our filaments number 1 and 2 are longer than the 
length considered by |Lee et al.|p014] ), while filament number 4 
has the same length. They identify the filaments based on their 
A 2 H^ emission, while our method is based on an enhancement 
of the column density due to the dust emission. For the three fil¬ 
aments in common, the masses given in [Lee et ar|p014| ) agree 
within the errors with the values we computed. Also the temper¬ 
atures are consistent with the temperatu re range we meas ured 
(see Table |^. We therefore agree with |Lee et al.] p014| ) that 
filament number 1 is thermally supercritical. 
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V. Roccatagliata et al.: Herschel observations of the Serpens Core. 




Fig. 2. Left: Temperature map of the Serpens Core in K. Right: Column density map (A^h) of the Serpens Core in The center of 

the image corresponds to the position (crj 2 ooo. ^moo) = (18^ 29^ 18% +00° 51' 08") and the size is 94' x 67'. In both panels the large boxes define 
the position of the whole Serpens Core. The nine filaments considered in our analysis are highlighted by numbered boxes. In filaments 1 and 7 the 
small squared boxes highlight the positions of the cores. Box 0 represents the control field. 


Table 2. Position, size, total mass, mass per length, critical mass per length, and temperature range derived for the filaments highlighted in 
Figure]^ For comparison we report the total mass and temperature range of the Serpens Core. 


Observations of the Serpens Core 


Filament 

Coordinates^^^ 

Size^^^ 

Angle*^^^ 

^dust+gas 

Ml 

^L,crit 

j(4) 

number 

0^12000 

^12000 

['] 

[pc] 

n 

[Mq] 

[Mopc-i] 

[Mo pc-i] 

[K] 

1 

Igh 30” 02** 

+01° 08'19" 

1.9X8.3 

0.2X1.0 

353 

32 

6:! 

32 

24 

12-17 

2 

IS** 30” 12* 

+01° 13'21" 

1.3X3.7 

0.2X0.4 

79 

13 

27 

13-19 

3 

IS** 30” 20* 

+01° 16'31" 

1.7X6.3 

0.2X0.8 

77 

d 

7 

27 

14-18 

4 

Igh 30” 13* 

+01° 19' 32" 

1.5X4.8 

0.2X0.6 

40 

5 

27 

15 - 17 

5 

18*’29” 57* 

+01° 24' 17" 

2.2x10.0 

0.3X1.2 

6 

7 

26 

15-17 

6 

18‘’29”43* 

+01° 20' 58" 

1.5X5.5 

0.2X0.7 

330 

13% 
23 % 

6 

25 

14-16 

7 

18** 29” 28* 

+01° 20' 07" 

1.7x10.4 

0.2X1.3 

319 

10 

24 

12-16 

8 

IS^ 28” 53* 

+01° 14' 06" 

2.2x17.4 

0.3X2.1 

68 

11 

26 

14-17 

9 

18** 28” 51* 

+00° 50' 09" 

2.9x14.3 

0.3X1.7 

327 

23% 

13 

24 

12-17 

Serpens core 

18** 29” 50* 

+01° 14' 14" 

42x42 

5x5 

36 

494 



18-30 


field 


The coordinates represent the center of the boxes. Size and orientation of the boxes. The angles are measured from North 
over East. Lowest and highest temperatures in the corresponding boxes. 


5.2. column density distributions 

The column density distribution of the entire field that covers the 
Serpens Core (as highlighted in Figure]^ is shown in Figurej^as 
a red histogram. The observations show a cut at densities lower 
than 10^^-^ cm“^, due to the sensitivity limit of the observations. 
The sensitivity limit is sufficient to well sample the peak in the 


observed column density distribution, which is more narrow than 
the peak in the simulations. Overplotted are the column density 
distributions from two simulation snapshots separated in time by 
0.4Myr. The distribution evolves relatively little during this 
period of time. 
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Fig. 3. Temperature map and column density map (A^h) of the high-resolution simulation of Run I of Dale et al. (20121. The boxes represent the 
regions which contain filaments and cores. Box 0: control field; Boxes 1 & 2: cores; Boxes 3,4 & 5: filaments. The small protocluster is contained 
in Box 1. 


Table 3. Mass and temperature range derived for the simulated fila¬ 
ments highlighted in Figurej^ For comparison we report the total mass 
and temperature range of the entire region that appears in Figure 


Simulation 

Region 

Size 

M 

T 

# 

pcxpc 

[Mo] 

[K] 

3 

0.3X0.7 

20 

19-23 

4 

0.3X1.2 

22 

21-28 

5 

0.4X1.3 

36 

22-27 

Entire Field 

5x5 

801 

23-50 


The slopes of the histograms of column densities have been com¬ 
pared between observations and simulations, as it is shown in 
Figures]^ to The histograms of column densities of the simu¬ 
lations and the PDFs have been constructed as in the analysis of 
the observations. The dashed-dotted lines in the left panel of Fig- 
ure|^represent the fit of a Gaussian distribution and a power-law. 
We find that the deviation from a log-normal distribution corre¬ 
sponds to Ay ~1.2-2.0 mag. Both column density distributions 
of the two snapshots of the simulation separated by 0.4 Myr in 
time, show a very broad peak resembling a log-normal distribu¬ 
tion, but with a pronounced tail at high column densities, rather 
different from that seen in the observed cloud. In particular, there 
is a considerable amount of low-density material that does not 
appear in the observed data (it is below the observational sensi¬ 
tivity limit), and the peak of both simulated distributions occurs 
at a higher column density. The quantity of high-column density 
material, between 10^^ and 10^^ cm“^ drops between the two 
simulation snapshots. The quantity of material which appears 
to be ‘missing’ between the two distributions is ~ lOM©, cor¬ 
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responding to an order of one percent of the total mass of gas 
in the simulation volume, and to about thirty percent of all gas 
denser than 10^^ cm“^. In the 0.4 Myr time interval between the 
two snapshots some of the high-density gas is consumed by star 
formation. The figure highlights the fact that, while a simulation 
snapshot and an observed field may appear subjectively similar, 
more careful analysis can reveal important differences between 
them. 

The distribution of the observed column densities of all fila¬ 
ments, including the cores, shows a peak at Ay ~ 0.6 mag, with 
a maximum value of Ay ~ 18 mag (see left panel in Figure [^. 
For Ay -1.5-14 mag, the distribution decreases as a power-law, 

following the relation dlog{N)ldlogiNH) ~ tr, where a -1.3 

(red histogram in the left panel of Figure [^. In the simulations 
the peak is found at a higher value i.e. Ay ~ 3 mag. The slope 

computed between Ay - 6 -12 mag is a -1.6, i.e. steeper than 

the observed value. 

The core emission has been identified as an enhancement in the 
column density map, as well a point source emission at least at 
one Herschel wavelength. In filament number 1 in the Serpens 
Core, we found two cores (both included in the same box) and in 
filament number 7 one core. Since all the cores do show a point 
source counterpart at 70 jam, this means that we are selecting the 
more evolved protostellar cores, and not the cores at an earlier 
evolutionary stage (when they do not have yet 70yum emission). 
In the simulations two cores are identified. The histogram of 
the filaments without cores, obtained by drawing a box around 
the cores and excluding them, is computed for both observations 
and simulations (see left panel of Figure |^. In the observations 
the maximum column density reaches Ay - 9 mag, and its peak 
is at Ay ~ 0.7 mag. While the observations show a peak with 
a power-law tail with a slope a ~ -1.5, the simulations show 
only a peak at Ay ~ 3 mag. The slope in the observed filaments 
becomes steeper excluding the cores. 

















































V. Roccatagliata et al.: Herschel observations of the Serpens Core. 


Av 

0.05 0.5 5.0 50 500 




- 2-10 1 2 

ri=ln(N /<N >) 

H H 


Fig. 4. The left panel presents the column density distributions and the right panel the normalised PDFs of the Serpens Core (in red) and 
two snapshots of Run I of the simulations of [Dale et al.| ( |2012] ) that differ by 0.4 Myr in time (in blue and green). The snapshot at tsc is shown 
in Figurej^ The dashed-dotted red lines in the left panel represent the fit of a Gaussian distribution and a power-law. In the right panel the 
dashed-dotted lines represent the best Gaussian fit of the observations (in red), Run I at tsc (in blue) and Run I at tsc+OA Myr (in green). 


We also did the analysis of the filaments with and without cores 
excluding the filament number 1, that is the only supercritical 
filament in the network of filaments of the Serpens Core. In 
Figures |^and [^the histograms without filament number I are 
drawn in green. In both cases the histogram still shows a power 
law tail, but with a steeper slope. In the case of the filaments with 
cores a steeper slope was expected, since by excluding filament 
1, we are not considering the two cores that are present in this 
filament. It is even more interesting to notice that even without 
the only supercritical filament, the observed filaments still show 
a power law tail. 

A histogram showing the cores alone has been constructed ex¬ 
cluding the filaments in both observations and simulations (see 
left panel of Figure [7]). In the case of the simulation the maxi¬ 
mum column density reaches a value of Ay ~ 12 mag, and the 

slope decreases to an almost flat slope of a -0.9. In this case 

of the cores alone there are not many pixels which sample the 
cores and the errors are much larger than in the histograms of 
the column densities of the filaments with and without cores. 


5.3. Probability density functions (PDFs) of the filaments. 

The PDFs of the same regions described in Section |5.2| have 
been computed for the observations and simulations and they 
are shown in the right panels of Figures]^ to |7] 

A quantitative comparison between the PDFs of the observa¬ 
tions and the simulations has been carried out using the standard 
non-parametric Kolmogorov-Smirnov statistic test. The com¬ 
puted probability that the observed and simulated distributions 
are drawn from the same distribution is give in Table We no¬ 
tice that the lowest probability is found in the case of filaments 
with cores. The value is, however, not sufficiently low to con¬ 


clude that simulations and observations are from different parent 
distributions. 

We perform a Gaussian fit of each PDF in order to characterise 
the log-normal function: 


Prjdrj = 


1 




exp 


I- 


(j] - nf 

20 -^ 


]di] 


(3) 


where cr^^ is the dispersion, and p is the mean logarithmic column 
density. The widths, cr^, resulting from the fitting are reported in 
Table]^ We notice that the values of cr^^ in the observations of the 
Serpens Core are always lower than the values in the simulations. 
In the case of the cores alone, the values of cr^ are not meaningful 
since both, the observed and simulated PDFs do not show any 
peak. 

As for the column density distributions, also the PDFs of the 
filaments alone in the simulations show only a log-normal dis¬ 
tribution while the observed filaments do show a power-law tail 
(see right and left panel of Figure |^. This is in agreement with 
the finding that the filaments themselves in the simulations are 
not fragmenting, but are instead serving as conduits of gas into 
the main cluster. The simulated cores at higher column densities 
show instead the power-law tail, while the observed cores have 
an almost flat distribution (see Figure [7]). In the PDFs the log¬ 
normal component always dominates compared to the power- 
law tail and the simulations show always a higher FWHM than 
the observations, apart from the case of the cores alone. Since 
the FWHM of the PDFs is a measure of the turbulent dispersion 
velocity (e.g. [Federrath & Kless^|20I3| ), this implies that the 
turbulent velocities in the observed filaments are always lower 
than those in the simulation. 

A power-law tail in the column density distribution has been in¬ 
terpreted as a transition between a turbulence-dominated gas and 
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Fig. 5. The left panel represents the histograms of the column density distributions while the right panel shows the normalised PDFs of 
the Filaments with cores. The dashed-dotted lines overplotted to the PDFs represent the best Gaussian fit of the data. Red lines represent the 
observations, the blue ones represent the simulations, and the green lines represent the observations without Filament number 1. The positions of 
the observed filaments and cores are highlighted in Figurej^ while the simulated filaments and cores are outlined in Figurej^ 
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Fig. 6. The left panel represents the histograms of the column density distributions while the right panel shows the normalised PDFs of the 
Filaments without cores. The dashed-dotted lines overplotted to the PDFs represent the best Gaussian fit of the data. The positions of the observed 
filaments are highlighted in Figurej^ while the simulated filaments are outlined in Figurej^ The different colors of the lines are as in Figure]^ 


a self-gravitating gas (e.g. |Klessen||2000t [Federrath & Kless^ 
|2013| ). Studies of the probability distribution functions, com¬ 
puted from the column density maps, have been used to constrain 
star formation processes. Using the column density maps from 
Herschel, [Schneider et al.| ( |2015| ) compared the PDFs of low- 
mass and high-mass star-forming regions, finding a common 
“threshold” value of Ay ~4-5 mag (corresponding to Ah =3.8- 


4.7x10^^ cm“^) between turbulence-dominated gas and self- 
gravitating gas. This “threshold” is higher than the value of 
Ay ~1.2-2.0 mag that we find in the column density distribu¬ 
tion of the Serpens Core (see Figure]^. The column density 
PDFs presented suggest that the simulated filaments are likely 
turbulence-dominated, while the observed filaments have the 
characteristic power-law tails indicating that they are most likely 
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V. Roccatagliata et al.: Herschel observations of the Serpens Core. 
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Fig. 7. The left panel represents the histograms of the column density distributions while the right panel shows the normalised PDFs of the cores 
alone. The positions of the observed cores are highlighted in Figure]^ while the simulated cores are outlined in Figurej^ The different colors of 
the lines are as in Figurej^ 


Table 4. Standard deviation cTj^ (defined in Equation ^ of the log¬ 
normal distributions overplotted in Figures|^[^|^&[^ Tne last column 
shows the probability that the simulated and observational data (the ob¬ 
served filaments refer to numbers 1 to 8) are drawn from the same dis¬ 
tribution using the K-S test. 



Serpens Core 

Run I @ tsc 

Pk-s 

Entire Field 

0.27 

0.78^ 

0.54 

Filaments with cores 

0.28'’ (0.25^^) 

0.36 

0.03 

Filaments without cores 

0.24'’ (0.25^^) 

0.32 

0.19 

Control field 

0.13 

0.38 

0.36 


(Tr^ of Run I at tsc+0.4 Myr is 0.68. Filaments numbers 

1-8. Filament numbers 2-8 (excluding Filament number 1). 


in the self-gravitating regime and that they will probably even¬ 
tually fragment into cores. A more stringent criterion for grav¬ 
itational instability which indicates that filaments are already in 
the process of sub-fragmentation (e.g. | Andre et al.|2010| 2Q14| ) 
is whether their masses per unit length are supercritical. Recent 
observations by |Lee et~ar] ( |2014| ) found that filaments FSl, FS2 
and FS3 (which correspond to our filament number 1) do indeed 
have supercritical values of M/L, while their filaments FCl and 
FN1 (our filaments 2 and 4, respectively) have M/L values close 
to the critical threshold (Table |^, suggesting that they will begin 
sub-fragmentation very soon. 


6. Summary and Conclusions 

In this paper we analyze continuum Herschel observations of the 
Serpens Core and compare them to a numerical simulation from 
the suite of calculations performed by [Dale et al.| ( |2Q12| ). Since 
the Serpens Core does not harbor any high-mass star, we choose 
an early state of a particular simulation of [Dale et al.| ( |2012) be¬ 


fore any 0-stars have formed and therefore without any feed¬ 
back. It is possible that this region will form high-mass stars 
in the future, since it has a substantial reservoir of gas available 
and its maximum column density is only a factor of a few lower 
than the canonical 1 g cm“^ threshold for massive star formation 
proposed by Krumholz & McK^ ( [2008| ). In the chosen simula¬ 
tion, the action of the filaments is mainly to deliver more gas into 
the cluster which forms at the hub where they join, without frag¬ 
menting. We selected a sub-region of Run I with the most simi¬ 
lar total mass, spatial scale and morphology to the Serpens Core. 
We then compared the properties, such as temperature and sur¬ 
face density in regions identified by the presence of a filamentary 
structure or a core. Note that the simulation was not designed to 
produce an object resembling Serpens (or any other star-forming 
region). The resem blance is serendipitous and owes to the aim of 
[Dale et aL] ( |2012| ) to cover a realistic parameter-space with gen¬ 
eral simulations. The picture which emerges is rather complex. 
With the additional time dimension available in the simulations, 
we saw that, while the simulations and observations may appear 
similar in morphology at a given time, they may differ in other 
respects, such as in the form of the column density distribution. 
Advancing the simulation in time resulted in the morphological 
disparity increasing, but the two datasets becoming more simi¬ 
lar from the point of view of their column density distributions. 
This indicates that quantitatively comparing simulations and ob¬ 
servations is more difficult than it may firstly appear. We defer 
a more detailed analysis to a later paper in which we construct 
synthetic observations of the simulations in order to make such 
a comparison more robust. 

The main results and conclusions of our study can be summa¬ 
rized as follows: 


1. Herschel continuum observations between 70 and 500 pm re¬ 
veal a radial network of filaments converging into the central 
cluster of the Serpens Core: this low-mass star forming re- 
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gion has been used to test the state-of-the-art of simulations 
ofGMCs. 

2. In the filaments the dust temperatures are found to lie in the 
range between 12 K and 19K in the observations and few 
Kelvin lower than the gas temperatures in the simulations: 
this might reflect a difference in temperature between gas 
and dust in the region. 

3. The peak in the observed column density distribution of the 
entire field of the Serpens Core is more narrow than the peak 
in the simulations. Both distribution show a power-law tail 
at higher column densities but with different slopes. This can 
be explained by a different evolutionary stage of the simula¬ 
tion in comparison to the real observations. 

4. The peaks in the column density distributions of the fila¬ 
ments are about 10^'^ cm“^ higher in the simulation than 
those in the observations. 

5. From the measured widths, cTr^, of the PDFs, we found that 
the turbulent velocities in the observed filaments are likely 
lower than those in the simulations. 

6. The column density distribution and PDFs of simulated fil¬ 
aments shows only a log-normal distribution, while the ob¬ 
served filaments show a power-law tail. The power-law tail 
becomes evident in the simulation in the column density dis¬ 
tribution of the cores, while in the observations this shows 
a fiat distribution. This means that while the simulated fila¬ 
ments are likely turbulence-dominated regions, the observed 
filaments are probably turbulence-dominated self-gravitating 
structures which will fragment into cores. 

In the simulation, the main function of the filaments is to 
feed gas into the central cluster which is where most of the 
star formation takes place. Most of the gas in the filaments 
falls into the cluster before fragmenting. It is therefore pos¬ 
sible that, despite the physical and morphological similari¬ 
ties between the chosen simulation timesteps and the Ser¬ 
pens Core, the filaments in the real and simulated clouds are 
behaving in very different ways. Whether this is due to ne¬ 
glected physics in the simulations (e.g. magnetic fields), to 
unknown differences between the early phases of the Ser¬ 
pens Core and the initial conditions of the calculations, or 
just to the extreme difficulty of comparing two very complex 
dynamical systems is not clear. The simulated and observed 
filaments are subjectively similar-looking, but behave differ¬ 
ently, the former being turbulence-dominated and the latter 
being self-gravitating. 
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Appendix A: The Control field 

We consider in both observations and simulations a control re¬ 
gion where neither filamentary structures nor cores are present 
(box 0 in Figures[^and[^. Figure [AT] shows the histograms of 
the column density distributions on the left, and the PDFs on the 
right. Since only few bins are covered in the column density dis¬ 
tribution, these plots are shown in the appendix. However, the 
small error-bars in each bin highlight that the column density 
distributions show only a peak without any additional structure. 
This peak at Ay ~ 0.5 mag is very narrow in the observations, 
while in the simulations the peak is much broader. 
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Fig. A.l. The left panel represents the histograms of the column density distributions of the control field, while the right panel the PDFs. The 
dashed lines overplotted to the PDFs represent the best Gaussian fit of the data. Red lines represent the observations, while the blue ones represent 
the simulations. 



Fig. B.l. Filament number 1. Starting from the left side, the panels show HerschelfPACS 10 jum, 160 jim, and SPIRE 250 yum, 350 yum, and 
500yum, respectively. The unit of all the images is Jy/arcsec^. The last two panels show the column density map (Nu) in [x 10^^^cm"^] and the 
temperature map in K. The box identifies the extension of the filament and its size is listed in Table The smaller box identifies the position of 
the cores among the filament. 



Appendix B: Filamentary structures in the Serpens 
Core. 

Figure [BT shows the enlargements on Filament number 1: start¬ 
ing from the left side, the first five panels show the Herschel 
PACS and SPIRE mosaics, while the last two panels show the 
column density and the temperature map. An additional box 
identifies the position of a core structure among the filament. 

In the online material, Figures show the enlargements on 
each filament as in Figure [BG] 
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Fig. 8. Filament number 1. Starting from the left side, the panels show HerschellVAC^ 70 pm, 160 pm, and SPIRE 250 pm, 350 pm, and 500 pm, 
respectively. The unit of all the images is Jy/arcsec^. The last two panels show the column density map (Ah) in [xl0^^^cm“^] and the temperature 
map in K. The box identifies the extension of the filament and its size is listed in Table The smaller box identifies the position of the cores 
among the filament. 
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Fig. 9. Filament number 2. The panels and the units are as in Figurej^ 




Fig. 10. Filament number 3. The panels and the units are as in Figure 
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Fig. 11. Filament number 4. The panels and the units are as in 
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Fig. 14. Filament number 7. The panels and the units are as in Figure The smaller box identifies the position of the core among the filament. 
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Fig. 15. Filament number 8. The panels and the units are as in Figure 
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Fig. 16. Filament number 9. The panels and the units are as in Figure 







































